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Abstract 

We use the SSTL (Singwi, Sjolander, Tosi, Land) approximation to investigate the 
short-range correlations in a one dimensional electron gas, for the first time. Although 
SSTL is introduced to better satisfy the compressibility sum rule in three dimensions, 
the widely used STLS (Singwi, Tosi, Land, Sjolander) approximation turns out to be 
more successful in the case of the one dimensional electron gas. 
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1 Introduction 



The advances in fabrication technologies have made it possible to experimentally fabricate 
one dimensional electronic structures This has naturally resulted in an increasing in- 
terest in theoretical investigations of such structures [H-flSl. The resulting one dimensional 
electron gas has also been subject to computational investigation [O . 



The extensive investigations in the three dimensional electron gas have clearly shown 
the importance of the short range correlations at lower densities in determining the pair 



distribution function g{r) at small r|T5 



In the high density limit, the long range correlations can be described well by the random 
phase approximation (RPA). One way to include the short range correlation effects beyond 
RPA is to use the powerful approach developed by STLS|ll6[|. 

It is well known that STLS approach suffers from a compressibility inconsistency. The 
compressibility calculated from the small q limit of the dielectric function e{q) does not 
agree with that calculated from the ground state energy. There have been several attempts 
to overcome this inconsistency [p!7|, [18|. SSTL|]l^ included the effect of screening on the 



effective interaction potential as explained in the next section. 

In this work, we aim to check the performance of the SSTL approach and compare it 
with the STLS in a one dimensional electron gas. To the best of our knowledge, this is the 
first application of the SSTL approach in lower dimensions. 

The organization of the paper is as follows; the STLS and SSTL formalisms are given 
in section 2. The results and discussion concentrating on the comparison between STLS 
and SSTL performance on the compressibility issue are presented in section 3. 

2 Formalism 

In the mean field approximation, a key component of the electron gas is the usual static 
structure factor S(q) which is related to density-density response function x(g, c<j) through 
the fluctuation-dissipation theorem as 

1 f°° 

Siq) = / duj Imx{q,uj). (1) 
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The density-density response function cu) is defined as 

1 - ^e//(g)xo(g,^) 

where Xoil:^) is the zero-temperature susceptibility of a noninteracting electron gas and 
is given in one dimension by 

Xo (<?, t^) = Xoi (<?, tu) + iXo2 (g, c^) , (3) 

with 



in 

Xoi(g,^) = 72 — In 
n Txq 



(4) 
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Ti q 



UJ- < UJ < UJ< 



otherwise, 



where uj± 



are the boundaries for particle-hole excitations. The Fermi 



2m* m* 

wave vector kp is related to the one dimensional electron density n via kp = /2. The 
dimensionless electron density parameter is defined as = / {Akpa*^), where a% is the 
effective Bohr radius. 

Veff{q) in Eq. (2) is the self-consistent effective potential related to the static structure 
factor through 



1 POO rjh 

Veffiq) = v{q) + ^ [Sik) - 1] {(g + k)v{q + k) + (q ~ k)v{q - k)} , (5) 

where v{q) is the Fourier transform of the Coulomb interaction between two electrons in 
the lowest subband in the harmonic confinement in one dimension and is given by 

v{q) = — Fiq), (6) 

^0 

where is the background dielectric constant, q is the wave vector along the wire and F{q) 
is the form factor which takes into account the finite thickness of the wire. The form factor 
reads 

(7) 

here Kq{x) is the zeroth order modified Bessel function of the second kind and b is the 
characteristic length of the harmonic potential. Indeed, it is a measure of the effective 
radius of the quantum wire[|T^ . 



In the STLS approximation. 



1 - v{q)[l ~ G{q)]xo{q,uj) 
where the local field correction G{q) is 

The set of Eqs. (1), (8) and (9) have to be solved self-consistently for G{q), and S{q) 
within the STLS approximation. In RPA, G{q) = 0. 

The SSTL approximation is different from the STLS approximation in that the potential 
under the integral sign in Eq. (5) is screened by the static dielectric function e{q) which is 
given by. 
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1 + G'(g)t;(g)xo(g) 

This is originally done to better satisfy the compressibility sum rule in three dimensional 
electron gas. 

The ground state energy per particle in one dimensional electron gas may be written 
as[gg] 

where Ekin is the kinetic energy per particle, and simply e^m = tt^/ (48r^) in units of effective 
Rydberg Ry*. Sex is the exchange energy per particle 

\ rkp rkp-q 

Sex = -7n~ dk v{k), (12) 

/TT^ra JO J~{kp+q) 

and Scot is the correlation energy per particle 

1 



dr', / dqv{q) [5(^,0-1]. (13) 
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The compressibihty K may be calculated either by using the ground state energy per 
particle Eg 



^-"^Isl-a (14) 

or by using the g ^ limit of the static dielectric function e{q) 

lim e{q) = 1 + v{q) K. (15) 

The compressibility sum rule is that the compressibilities calculated by using Eq. (14) 
and Eq. (15) are the same. 



3 Results and Discussion 

In Fig. 1 the local field correction G{q) calculated self-consistently using STLS and SSTL 
approximations are shown for fixed b = 5a^ and = 5. The most striking difference 
between the two curves is their large q limits. As STLS G{q) is approaching 1, SSTL G{q) 
approaches a limit bigger than 1. It may be noted that the relation[|r^ G{oo) = 1 — g{0) 
is satisfied in our calculations. This, of course, is going to lead to different results when 
applied to physical properties of the one dimensional electron gas. 

The different large q behaviour of G{q) is expected to lead to different small r behaviour 
of the pair correlation function g{r). This is shown in Fig. 2. As radius of the wire gets 
smaller the difference between STLS and SSTL results become larger, as may be seen in 
Fig. 3. As the density becomes smaller (i.e., becomes larger) the SSTL g{r) becomes 
negative for small r. This limit is where the correlation effects become important. 
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The SSTL s{q)~^ is given in Fig. 4. This is rather similar to STLS or RPA e{q)~^ with 
a discontinuity in the derivative aX q = 2kp. 

The ground state energy per particle for our system is shown in Fig. 5. The STLS and 
SSTL curves are rather similar. 

The compressibility calculated by using the ground state energy in three different ap- 
proaches is presented in Fig. 6 as a function of r^. The similarity in energies is also 
reflected in these curves. The negative compressibility values at higher show that the 
system becomes unstable. 

The compressibility calculated by using g ^ limit of the static dielectric function is 
shown in Fig. 7. Here, the STLS and SSTL compressibilities calculated by two routes 
are as different as in RPA. This surprising result is a natural consequence of the small q 
behaviour of the local field correction which is very small for SSTL at small q values. 

The compressibility in a one dimensional system is studied previously by Gold and 
Calmels|0] using a three-sum-rule approach for the local field correction. The compressibil- 
ity sum rule is better satisfied within this variant of STLS in one dimension as can be seen 
from their figure 4. We find the same result by the present full STLS calculation. 

We can conclude that the short-range correlations are calculated using the SSTL ap- 
proximation in a one dimensional electron gas, for the first time. The performance of the 
SSTL approximation is compared with the more widely used STLS approximation. It is 
shown that SSTL compares reasonably well with STLS when the pair correlation function, 
dielectric function and the ground state energy per particle are considered but fails totally 
when the compressibility sum rule is checked. 
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Figure 1: The local field correction in SSTL and STLS approximations for wire radius 
b — 5a% and = 5. 




Figure 2: The pair correlation function in SSTL and STLS approximations for wire radius 
b — 5a^ and = 5. 
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Figure 3: The pair correlation function in SSTL and STLS approximations for wire radius 
and Tg 



b — 2a% and r, = 5. 




Figure 4: The inverse dielectric function in different approximations for wire radius b — 5a^ 
and r, = 5. 
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Figure 5: The ground state energy per particle in different approximations for wire radius 
b = 2a%. 
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Figure 6: The compressibihty calculated by using the ground state energy in three different 
approaches for wire radius b = 2a^. Hereafter, Kq is the free electron compressibility. 
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Figure 7: The compressibility calculated by using the lirn e{q, 0) in three different 
proaches for wire radius b = 2a%. 
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